Microscopic approach to orientational order of domain walls 
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We develop a fully microscopic, statistical mechanics approach to study phase transitions in Ising systems 
with competing interactions at different scales. Our aim is to consider orientational and positional order param- 
eters in a unified framework. In this work we consider two dimensional stripe forming systems, where nematic, 
smectic and crystal phases are possible. We introduce a nematic order parameter in a lattice, which measures 
orientational order of interfaces. We develop a mean field approach which leads to a free energy which is a 
function of both the magnetization (density) and the orientational (nematic) order parameters. Self-consistent 
equations for the order parameters are obtained and the solutions are described for a particular system, the Dipo- 
lar Frustrated Ising Ferromagnet. We show that this system has an Ising-nematic phase at low temperatures in 
the square lattice, where positional order (staggered magnetization) is zero. At lower temperatures a crystal- 
stripe phase may appear. In the continuum limit the present approach connects to a Ginsburg-Landau theory, 
which has an isotropic-nematic phase transition with breaking of a continuous symmetry. 

PACS numbers: 64.60.De,75.70.Ak,75.30.Kz,75.70.Kw 



I. INTRODUCTION 

Systems with competing interactions are a rule in nature. 
Competing tendencies are usually responsible for the com- 
plex behavior of natural systems, leading to slow dynam- 
ics, metastability and complex free energy landscapes, like 
observed in spin glasses and many other frustrated systems. 
Competing interactions at different scales may give rise to 
complex phases and patterns, like stripes, lamellae, bubbles 
and others^. Examples range from solid state systems, like ul- 
trathin ferromagnetic films ^ and strongly correlated electron 
to soft matter systems like Langmuir monolay- 



systemsi^. 
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ersS, block copolymers^, colloids and soft core systems^ 
Besides the intrinsic interest raised by the complexity of the 
phase behavior in these systems, their detailed knowledge 
may be critical for understanding basic phenomena as high 
temperature superconductivity, and also for technological ap- 
plications like soft matter templates for nanoscale systems and 
future spintronic devices. 

Systems with competing interactions at different scales, e.g. 
a short range attraction and a long range repulsion, as present 
in magnetic films and low dimensional electronic systems, 
have a tendency to form microphases, i.e. to phase separate 
at mesoscopic scales, due to the frustration usually present 
as a consequence of competing effects. These microphases 
can show at least two types of ordering: positional ordering 
of the microscopic degrees of freedom and also orientational 
order of interfaces or domain walls, due to the presence of in- 
terfaces at mesoscopic scales. The presence of orientational 
order allows for an analogy with liquid-crystal phenomenol- 
ogy, where nematic, smectic and crystal phases can in princi- 
ple be characterized. Nevertheless, a complete characteriza- 
tion of the phases present in these systems is rare, and usually 
only positional order parameters are computed and the cor- 
responding phase diagrams are well known — ~. In stripe 



forming systems, like the ones studied in the previous refer- 
ences, these mean field approaches give access only to mod- 
ulated crystalline phases, or stripe-crystals, where both posi- 
tional and orientational long range order are simultaneously 
present. In fact, the morphologies of domains observed in real 
films, like the magnetization patterns in ultrathin ferromag- 
netic film o 2 ' 3 ' 17 for example, give clear indications that ori- 
entational order sets in well before any positional order can 
appear. It is necessary to go beyond the usual mean field 
approaches in order to describe the complete ordering pro- 
cess. To our knowledge, there have been very few attempts 
to compute full phase diagrams, with simultaneous and inde- 
pendent consideration of orientational and positional orders. 
One of these few works is the phase diagram for ultrathin 
ferromagnetic films with square lattice anisotropy anticipated 
by Abanov and coworkers^. These authors have analyzed 
the phases of an Heisenberg spin system in the square lattice 
with competing exchange and dipolar interactions, by means 
of a mixture of microscopic calculation with phenomenolog- 
ical assumptions. Their analysis gave a very rich phase dia- 
gram, with isotropic, Ising-nematic and smectic phases. Phe- 
nomenological approaches are common ground for the anal- 
ysis of orientational order. Starting from an elastic energy 
which assumes a crystalline ground state, specific fluctuations 
of the ground state can be studied perturbatively. A complete 
analysis of two dimensional smectic elasticity has been done 
thirty years ago by Toner and Nelson 



Another interesting approach to study the phase transitions 
in systems with isotropic competing interactions was due to 
Brazovskii — . Analyzing a generic Ginsburg-Landau model 
whose main characteristic is the presence of a minimum in 
the spectrum of Gaussian fluctuations at a non-zero wave vec- 
tor, he showed, in the context of a self-consistent field ap- 
proximation, that the model has a first order transition to a 
modulated phase, induced by fluctuations. While concep- 
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tually important, this result suffers from the same problems 
as other mean field approaches, namely, it describes phases 
where both orientational (modulated) and positional order sets 
in simultaneously and has very strong transversal fluctuations 
in two dimensions 2 ^. In a recent work we have gone beyond 
the Brazovskii approximation, by considering other terms in 
the Ginsburg-Landau free energy of two dimensional systems, 
which are dictated by the symmetry of the interaction s 22 ' 23 . 
Interestingly, pure symmetry considerations lead to terms in 
the free energy that encode orientational order parameters . We 
were able to find a nematic phase, at temperatures above the 
original Brazovskii modulated phase. By analyzing fluctua- 
tions of the mean field nematic solution we further found that 
the isotropic-nematic phase transition in the continuum sys- 
tem with isotropic interactions is of the Kosterlitz-Thouless 
type^ 3 -. These works set the stage for a complete Renormal- 
ization Group treatment of this kind of systems. Neverthe- 
less, universality in two dimensions is a delicate matter, and it 
would be important to get contact with the phenomenological 
theory from more microscopic, interaction specific, models. 

In this work we give a step further in the description of 
nematic order in systems with competing interactions, and 
develop a fully microscopic statistical mechanical approach 
which allows to consider positional as well as orientational or- 
der parameters on equal grounds. We study Ising systems with 
arbitrary competing interactions on a square lattice. We intro- 
duce a nematic order parameter on the lattice, by analogy with 
the order parameter studied in the framework of the Ginsburg - 
Landau model 2 ^. Then we develop a mean field approach 
which includes both positional (magnetization) and orienta- 
tional (nematic) order parameters in the free energy. Due to 
the two-body nature of the orientational order parameter it is 
necessary to go beyond one site approximations in order to 
compute the mean field free energy. A set of self-consistent 
equations for the order parameters, equivalent to the Bragg - 
Williams approximation for the Ising model, is obtained and 
the solution for a model system presented. We computed the 
phase diagram of the Dipolar Frustrated Ising Ferromagnet 
(DFIF), a well known and very studied model for ultrathin 
magnetic films with perpendicular anisotropy^i^i We show 
that this model has an Ising-nematic phase in the square lat- 
tice, where only orientational order is present. At lower tem- 
peratures a further transition to a modulated phase with po- 
sitional order, the stripe phase, is possible. In the continuum 
limit our approach leads to the Ginsburg-Landau model dis- 
cussed before 22 . It is also shown that, when considering the 
nematic order parameter, the magnetic (spin-spin) susceptibil- 
ity diverges at a higher temperature compared with the usual 
mean field calculation, in which the orientational terms in the 
free energy are absent. 

The paper is organized as follows: in section HH we intro- 
duce the nematic and the Ising-nematic order parameters. In 
this section we also develop a mean field approach to deal 
with general competing interactions in the square lattice. In 
jfllll we apply our technique to the Dipolar Frustrated Ising 
Ferromagnet. Finally we discuss our results and conclusions 
in gIVj 




FIG. 1. Image of domain patterns in a Fe/Cu(001) ultrathin film, 
taken by Scanning Electron Microscopy with Polarization Analysis 
(SEMPA). Courtesy of Prof. Danilo Pescia, ETH, Zurich. 



II. MICROSCOPIC NEMATIC-LIKE ORDER 
PARAMETER 

Figure Q] shows a typical pattern of magnetic domains in an 
ultrathin ferromagnetic film of Fe/Cu(001) with strong per- 
pendicular anisotropy. In this system, the competing interac- 
tions are the short range exchange interaction and the long 
range dipolar interaction. It is clear that this pattern does 
not present positional order, but it is possible to distinguish 
some degree of orientational order in the stripe-like pattern. 
We want to define a suitable order parameter for this kind of 
order and quantify it. 

Domain walls are observed at the transitions between pos- 
itive and negative values of the perpendicular magnetization. 
In an appropriate mesoscopic scale, it is possible to define a 
continuum magnetization density $(:?). The gradient of this 
quantity 

N(x) = V<f>(Z) = (d x <f>,d y <f>) (1) 

naturally defines a director that quantifies the degree of orien- 
tation of domain walls. 

However, the kind of order we are looking for is insensitive 
to the vector orientation, as shown in figure |2] Therefore, in 
analogy with the nematic order parameter of liquid crystals, it 
is possible to define a local tensor order parameter: 

Q l3 (x) ee 0(f) Udj - \d 2 5^ 0(f), (2) 

where i, j = {x, y}. This symmetric and traceless tensor has, 
in two dimensions, only two independent elements, which es- 
sentially represent the mean orientation of domain walls and 
the strength of the orientational order. Similarly to the vector 
order parameter, a non zero value of Qij represents a break- 
down of the rotational symmetry 0(2). However, differently 
from a vector, this order parameter is invariant under rotations 
by 7r characterizing a nematic symmetry. 

A comprehensive description of the physical meaning of 
this order parameter can be found in ref. |28l Since the tensor 
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FIG. 2. (Color online) Schematic representation of a stripe phase 
with a dislocation. The arrows represent the director field, which is a 
local measure of orientation of domain walls. 



is symmetric, it can be diagonalized and in the principal axis 
base it has only one independent component, which can be 
written as: 

(Q) = J d 2 kk 2 cos (261) C(k), (3) 

where C(k) is the structure factor of the system and k — 
fc(cos#,sin#). This expression makes possible the com- 
putation of (Q) from experimental or numerical simulation 
data, since the structure factor is a usual observable. In this 
way, the orientational order parameter quantifies the degree of 
anisotropy of the domain pattern (i. e. the anisotropy in C(k)). 

In ref.— it was shown that this nematic-like order parameter 
emerges naturally from symmetry considerations in the Lan- 
dau expansion of the free energy of a system with isotropic 
competing interactions. At mean field level plus Gaussian 
fluctuations 23 , the phase transition turns out to be in the 
Kosterlitz-Thouless universality class. This seems to be a ro- 
bust result, and was originally obtained in a somewhat differ- 
ent context by Toner and Nelson^. 

Due to the ubiquity of orientational phases, intermediate 
between crystal and fluid ones, in a great variety of condensed 
matter systems, it is important to have a unified microscopic 
description. This would also allow us to make closer contact 
with experiments and simulations. With this motivation, and 
also trying to clarify the nature of the nematic-like interaction 
term in the Ginsburg-Landau free energy analyzed in refs . 22 ' 23 , 
in this paper we develop a general classic statistical mechanics 
approach to obtain the phase diagram of generic competitive 
systems with Ising variables. 

Consider the following Ising Hamiltonian on a square lat- 
tice: 

ii i^T./,,.s,.s\ Y^n.s, (4) 

i,3 i 



where {Si = ±l,i = 1 . . .N} are N Ising variables, the 
sum runs over all i, j pairs in the lattice and Bi is an external 
magnetic field. The interaction matrix Jy is assumed to have 
a ferromagnetic part, Jy = J < if i, j are nearest- neighbor 
sites (NNS) and an antiferromagnetic part Jjj > if i,j are 
not NNS. These are the only constraints in order to apply the 
method. 

Consider a discrete version of the nematic order parameter 
defined in coupled to a conjugate external nematic field 
through: 

hf (^A a S k A b S k - S -f(AS k ) 2 ^j (5) 

where k is the lattice index, a,b = x,y and A a are lattice 
derivatives. The field /i£ b has the same symmetries as the ten- 
sor order parameter, then we can choose the coordinate system 
along the principal axes and write the conjugate field as : 

( h k \ 
K b =[ (6) 
V o -h k J 

The coupling term of equation (O can be written as 

^^2h t SiK^Sj, (7) 

where the matrix K is given by 

( +1 if j = i±x 
Ku = \ (8) 
[ -1 if j = i±y 

x and y are unit vectors along the x and y axis of the square 
lattice. In this way, the Hamiltonian now reads: 

n = \ zZ ( J « - hiKi ^ 5 ^ - S BiSi (9) 

The local orientational order parameter is defined as : 

p ohi 

where Z = X){s} ex P { — /^} i s the canonical partition func- 
tion, with f3 = 1/ksT. This order parameter has only one 
component and not two, like the one of eq. d2}. If (Ni) is 
positive, the director points along the x direction wile if it is 
negative, the director mainly points in the y direction. These 
are the only possible directions of the director. For this rea- 
son, if (Ni) 7^ 0, the resulting anisotropic phase is called 
Ising-Nematics, since it breaks the rotational point group of 
the lattice and it is invariant under rotations by tt. Along the 
paper, for brevity, we generally use the term "nematic" to re- 
fer to this phase, however, whenever we deal with a square 
lattice model, "Ising-Nematic" should be understood. 
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The global order parameter can be written as 

Q = \Y, K ^SiSj). (ii) 

Equation ( fTTT i is completely analogous to the continuous ver- 
sion given by equation ||3}. The anisotropic matrix K for the 
Ising-Nematic phase, play the same role of k 2 cos 26 for the 
nematic one. Thus, in the same way that eq. (O measures 
the degree of anisotropy of the structure factor C(k), the or- 
der parameter dTTt describes the degree of anisotropy in the 
nearest-neighbor correlation functions between the x and y 
directions of the lattice. A slightly different form of this or- 



der parameter has been used before in simulations of stripe 
forming system s 26 ' 29,30 . In the present work, we attempt to 
compute it in a statistical mechanics framework. 

The technical problem of computing the order parameter re- 
duces to the calculation of nearest-neighbor correlation func- 
tions in a proper approximation. In this work, we will de- 
scribe a mean-field like approximation based on the use of a 
Hubbard-Stratonovitch (HS) transformation. 

Introducing a real variable on the lattice where i is the 
lattice index) by means of a HS transformation on the original 
Hamiltonian ©(see, e.g. ref. l3lt) and exactly summing up the 
Ising variables Si, we obtain an effective Hamiltonian given 
by 



H 



eff 



[{*}] = i E *"<- / <<*"< - \ E B ^ -r E lo s cosh ( E v.- ] 

ij i i \ : / 



(12) 



Using the general relation between correlations of the original 
discrete and continuous variables^: 



(13) 



for the cases of interest in this work, of an isotropic interaction 
matrix Jy, the order parameter reads: 



(14) 



A. Mean field approximation for the order parameter Q 

To compute Q given by equation ( TBI ), we begin by consid- 
ering the partition function 



Z = Af V$ e -^//( $ ) 



(15) 



where Af is a normalization constant and H e ff(Q) is given 
by eq. (TT2"i >. We want to introduce an order parameter with 
nematic symmetry. For this purpose we introduce a symmetric 
traceless tensor Q^' , where k is the lattice index and a,b = 
x, y and write the partition function as 



Z = Af' J P$ e -^.//(*) J VQ e -0£»3V»S) 



(16) 



where VQ = J\ k dQ k ' b . Note that, in the square lattice, 
the introduction of a fully symmetric traceless tensor of rank 
two is redundant, since the Ising-Nematic order parameter has 
only one independent component. However, we prefer to use 
this notation to stress that the method is general and can be 
also used to treat continuous systems with 0(2) symmetry, 



where the full tensor is needed to describe the phase transi- 
tion. 

The Gaussian integral in Q is simply a constant which can 
be absorbed in the normalization factor Af' . Shifting Qk — > 
Qk — A r fc($), keeping the measure VQ invariant we have, 

Z=Af' Jv$e-P H °ffW J VQ e-^THQi-Ntf t 

(17) 

where the symmetric traceless tensor Nk(&) will be defined 
later. 

At mean field level, Q£ b is given by the saddle-point equa- 
tion: 



dlnZ 



which reads 



dQl 

Qt = W> 



(18) 



(19) 



The average is given by: 



(jV£ 6 ) = 1 f V<5> N£ b (<I>) e -P( H -ii^)+T. k Tr(Q k -N k f 

z{Q) J 

(20) 

Choosing the coordinate system along the Qf^ principal 

axes 



and 



Qk b = Qk 



1 

-1 



1 

-1 



(21) 



(22) 



and defining ^2 k Nk = § Ylij Kij<&i<&j, the effective Hamil- 
tonian takes the form: 
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H [{*<, Qi}] = \Y1 ** J « $ i ~ \ E lo § cosh ( E J *i*i 

ij i \ J / 

+ 2 ^(/V 4 ) 2 + 2 ^ Q 2 - 4 ]T Q,/V l5 



(23) 



where the external magnetic field i?i has been set to zero. This 
is a non-quadratic Hamiltonian in the variables <!>;, in the pres- 
ence of a "mean field" orientational order parameter Q. The 
next step is to integrate the partition function in $>i considering 
quadratic fluctuations of around a saddle-point approxima- 
tion. 

Thus, the expansion of H up to quadratic order in the $ 
fluctuations reads: 



and 5® = $ - $ s 



In this approximation (5$i5$j) = jgiH")^ and then the 
mean field equation for Q, eq. dT4b . reduces to: 



(25) 



= i/($5P) + ^E ff «l*=* s , C 24 ) 



where $sp is the saddle-point solution H'(&sp) — 0, H'l 



In order to get the desired self-consistent equation for Q, 
it is necessary to find an explicit expression for iJ"($sp). 
After a lengthy but straightforward calculation, we obtain for 
the Hessian: 



m\SP 



■Jim -0^2 



JilJii 



QK hn + 



16 r 



SP 



K kl K k 



(26) 



Inserting this result in (125) a self-consistent equation for Q is 
obtained. It depends on the field $ SP and its square, which 
can be computed within the same mean field approximation. 



reads, in matrix notation: 



The previous result allows in principle the calculation of 
positional (<£>) and orientational (Q) order parameters. We 
expect that a purely orientational phase, with Q / and 
$ SP = 0, which will be identified with a nematic phase 
with orientational order of domain walls, but without posi- 
tional order, may be present for a class of systems with com- 
peting interactions. Considering this situation, the simplest 
self-consistent approximation for the nematic order parameter 

I 



16/3 



Tr 



K 



J - 2/3 J 2 - 2QK 



(27) 



This self-consistent equation for Q is the analog of the Curie- 
Weiss approximation for the magnetization in the Ising model. 
This is one of the main results of our paper. If equation (|27t 
has a non-trivial solution Q ^ 0, then the system exhibit an 
anisotropic but homogeneous phase with nematic symmetry. 
The presence, or not, of this phase depends of the detailed 
structure of the competing interactions, coded in the explicit 
form of the matrix J. 

Guided by the results of the continuum Ginsburg-Landau 
model 2 ^, we search for a critical point signaling an isotropic - 
nematic transition: Defining A = J — 2/3 J 2 and expanding 
the r.h.s. of ([27) in QTr(K/A) «lwe have 



16/3Q 



(28) 



If A is rotationally invariant (invariant under discrete rotations 

in the square lattice), then Tr (f ) = Tr (f ) 3 = 0. Then 
Q = is always a solution of the saddle point equations. If 



Q^O: 



( J 



(29) 



6 



and then, for Q <C 1: 



A. Nematic phase 



Q 



8/3 



1/2 



(30) 



This result implies a continuous, second order isotropic- 
nematic transition, at a critical temperature given by: 



(31) 



The existence of a solution of eq. (l3Tl l depends on the detailed 
structure of the matrix A. In the next section, we will show an 
explicit calculation on a model system in which this transition 
is present. 



III. THE DIPOLAR FRUSTRATED ISING 
FERROMAGNET 

The dipolar frustrated Ising ferromagnet (DFIF) is a sim- 
ple model for studying the thermodynamic phases of ultra- 
thin ferromagnetic films with strong perpendicular anisotropy 
2 ] 3 i 25 ~ 2Z. In the strong anisotropy limit, the anisotropic term 
of the classical dipole-dipole interactions is zero because the 
dipolar moments point perpendicular to the lattice plane, and 
the system degrees of freedom are well represented by Ising 
variables pointing perpendicular to the plane of the lattice. 
Furthermore, the strong uniaxial (perpendicular) anisotropy 
renders the energy symmetric with respect to rotations on the 
lattice. The Hamiltonian with competition between short- 
range exchange and long-range dipolar interactions can be 



wntten a 



,26. 



H 



J 
2 



E^ + fE 



g SiSj 



(32) 



<hj> 



The first sum runs over all pairs of nearest neighbor spins, and 
the second one over all pairs of spins of the lattice; is the 
distance, measured in lattice units, between sites i and j. The 
relevant parameter is the ratio between the exchange J > 
and the dipolar g > intensities. Then, we fix g = 1 without 
loosing generality. Note that the ferromagnetic, short range 
exchange interaction is frustrated by the long range, antifer- 
romagnetic, dipolar interaction. The possibility of an Ising 
nematic phase in these systems has been anticipated theoreti- 
cally by Abanov et al. — , and numerical evidence from Monte 
Carlo simulations has been reported inSH Several experimen- 
tal works have reported results which show domain patterns 
that could be identified with a nematic phased, but up to now 
its characterization and properties have not been discussed. 
Also, several theoretical, mainly numerical works have shown 
results of orientational order parameters'^, but no quantita- 
tive characterization or distinction between, e.g. stripe, smec- 
tic and nematic phases have been attempted so far. 



We have numerically solved equation (|3TT i for the DFIF 
model d32l , in the thermodynamic limit where the linear size 
of the system L — > oo. In reciprocal space, eq. (1311 reads: 



1 



d 2 k 



K(k) 



J(k)(l - 2p c J(k)) 



(33) 



where J(k) and K(k) are the Fourier transforms of the matri- 
ces Jij and Kij . T 
lattice is given by 



ces Jij and ify. The anisotropic function K(k) in the square 



K(k) = 2(cos k x a — cos k y a) 



(34) 



where a is the lattice spacing. The Fourier transform of the 
interaction function is given by 



J(k) = - J L nn (k) + g L dip (k). 



(35) 



where the the Fourier transform of the nearest-neighbor in- 
teraction L nn (k) = 2J(cosk x a + cos k y a) and the dipo- 
lar interaction for small k is approximated— by Ldi P (k) ss 
1 - ka+ (fca) 2 /4 

We have solved equation d33l in two different limits: keep- 
ing the full structure of the functions K(k) and J(k) where 
the lattice symmetry is preserved, and keeping only the lead- 
ing order terms in ka <C 1, where the effects of the lattice 
are suppressed (continuum approximation). The phase dia- 
gram in the T, J plane is shown in figure[3]for the two cases 
considered. Both cases give similar results in this mean field 
approach. The temperature scale depends on the value of the 
cutoff needed to get convergence of the integrals. This is a 
well known limitation of the HS transformation? 6 when the 
interaction matrix is not positive definite, as in this case. In 
any case, as the approximation is of a mean field nature, we 
only expect that the phase diagram be qualitatively correct. 

An isotropic-nematic transition without positional order, 
($} = 0, is obtained in both cases. It is important to stress that 
in the full lattice calculation, the low temperature phase is in 
fact an Ising-nematic, while in the continuum limit the phase 
is a nematic one. This last case corresponds to the result found 
in ref|H from a phenomenological Ginsburg-Landau model. 
A fundamental difference between both cases arise in the na- 
ture of the transitions. While in the continuum case the transi- 
tion is of the KT type, in the first case the universality class is 
probably Ising (see ref. Qj§for a discussion of different scenar- 
ios when lattice anisotropy is considered), but this remains to 
be proved for this model. The present microscopic approach 
for the nematic phase has the Ginsburg-Landau Hamiltonian 
studied in refs^2^ as the continuum limit. Note that the con- 
tinuum limit of equation d27l i corresponds to equation (12) of 
ref. SI 

Interestingly, orientational order develops at higher temper- 
atures for systems where the competing interaction is weak, 
i.e. for weak frustration. This corresponds to the experimental 
situation in ultrathin ferromagnetic films, where visual inspec- 
tion of the domains formed point to the existence of a nematic 
phase^ 
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FIG. 3. (Color online) Isotropic-nematic transition line in the dipolar 
frustrated Ising ferromagnet. 



In figure |4] we show the order parameter Q near the transi- 
tion, obtained solving eq. ( T30T > for the lattice case. It shows 
the second order nature of the phase transition. This is an ex- 
pected result since, in two dimensions, the nematic symmetry 
only allows even powers of the order parameter (nematic as 
well as Ising-Nematics) in the free energy^. This is different 
from three dimensions where the transition is of first order. 
However, as was already indicated, fluctuations could proba- 
bly drive the transition to a different class. 




FIG. 4. (Color online) Mean field nematic order parameter for J = 
2. 



B. Positional order 

At lower temperatures, besides the orientational order rep- 
resented by the nematic phase, positional order may also 
emerge. Within the present framework, the positional order 
parameter is given by the density field <f>i, which can be com- 
puted together with Q(T) from equations (|23T > and d26*l i. We 
have not pursued this in the present work, which is devoted 
to the description of the nematic phase. Instead, it is interest- 
ing to analyze the behavior of the two-point spin correlation 



function in the nematic phase, which gives information of the 
possible growth of positional order at high temperatures. 

The spin-spin correlation function, eq. ([13), in reciprocal 
space reads: 



G(k) 



1 



1 



J(k) - 2/3 J(fc) 2 - 2QK(k) J(k) 



(36) 



This equation must be solved self-consistently with eq. d27l >- 
If the term proportional to Q(T) in the denominator is dis- 
regarded, we end with the usual, mean field approach to the 
computation of modulated phases. In figure [5] we show the 
function G(k) for Q = 0, for a temperature slightly above 
the mean-field critical temperature, which in this case is given 
by T c = maxzJ{k). The high temperature profile of G(k) 
is isotropic and a phase transition to a modulated phase with 
characteristic wave vector fcrj ^ takes place through a 
breaking of a continuous rotational symmetry. This is the 
"Brazovskii's scenario" for the isotropic-stripes phase transi- 
tion in systems with nearly isotropic competing interactions 20 . 
Instead, upon inclusion of the nematic order parameter cor- 




FTG. 5. (Color online) Structure factor of the DFIF near the 
paramagnetic-stripes phase transition, without considering the ne- 
matic correction term, computed from eq. (136) for J=10. 



rection, the spectrum above the low temperature modulated 
phases changes in an essential way. In fact, because the ne- 
matic transition breaks rotational symmetry, the structure fac- 
tor is anisotropic in this phase, as shown in figure [6] again 
for a characteristic temperature just above the transition to a 
modulated phase with positional order (divergence of a stag- 
gered magnetic susceptibility). Note that the spectrum has a 
broad support in the plane of wave vectors, a clear indication 
that there is no positional order in the nematic phase. Nev- 
ertheless, it shows two well defined maxima at characteristic 
wave vectors ±fco> which show the 7r-rotation symmetry of the 
nematic phase. Furthermore, both peaks are on the x-axis, a 
consequence of the breaking of continuous symmetry already 
imposed by the square lattice. This is the structural signature 
of an Ising-nematic phase. 

Another interesting question regards the 
(in)commensurability of the characteristic wave vectors. 
In this solution, as expressed e.g. by equation (l3o*l l. it is clear 
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FIG. 6. (Color online) Structure factor of the DFIF in the nematic 
phase, near the nematic-stripes phase transition for J=10. 




FIG. 7. (Color online) Structure factor of the BNNNI model near the 
paramagnetic-stripes phase transition, computed from eq. d36t and 
(37). 6 = 1.8 and /3 = 1/13. 



that the wave vector which maximizes the structure factor 
in the nematic phase depends continuously on temperature. 
Then, at least in the nematic phase, the characteristic wave 
vectors are incommensurate with the lattice. 

A very different situation arises in short-ranged inter- 
action models like anisotropic next-nearest-neighbor Ising 
(ANNNI) 39 or the biaxially next-nearest-neighbor Ising model 
(BNNNI)2P.. This class of models could be similarly treated 
within our formalism. For instance, the Fourier transform of 
the interaction matrix Jy for the BNNNI model reads, 

J(k) = 26 (cos(k x a) + cos(k y a)) 

- 2 (cos(2k x a) + cos(2k y a)) (37) 

where a is the lattice constant and 5 measures the competition 
between first and second neighbors interactions. In figure (Jvj 
we show the structure factor of this model, computed form eq. 
(f36l > in the high temperature phase. We observe the presence 
of four peaks obeying the lattice symmetry. The peaks weights 
grow as the temperature is lowered , signaling the onset of po- 
sitional order. However, the position of the peaks in recipro- 
cal space are given by cos(k x a) = cos(k y a) = 6/2 which 
do not depend on temperature. This fact allows a possible 
incommensurate as well as a commensurate positional order 
at low temperatures. Concerning the main focus of this pa- 
per (the nematic phase), we have computed the self-consistent 
equation PIT ) for the nematic critical temperature and we have 
found no solution for this case. This is consistent with the 
general belief that it is necessary to have a macroscopically 
degenerate number of ordering wave vectors to produce a pure 
orientational ordered phase (compare fig. (|5) with fig. ([7}). 

IV. CONCLUSIONS 

We have developed an approach to study phase transitions 
in two dimensional Ising systems with competing interactions 
which may show orientational and positional order. A nematic 
order parameter in the square lattice has been defined, suitable 
to quantify the degree of orientational order of interfaces, use- 
ful in systems which show microphase separation. A mean 



field approach has been developed which leads to a set of self- 
consistent equations for orientational and positional order pa- 
rameters, similar in spirit to the Bragg-Williams approxima- 
tion. The approach is very general, and can be applied in prin- 
ciple to any Ising system with competing interactions. Ne- 
matic, smectic and stripe-crystal phases can be studied within 
this framework. We solved the self-consistent equations for 
the Dipolar Frustrated Ising Ferromagnet. This model for 
ultrathin ferromagnetic films with perpendicular anisotropy 
is known to have a striped low temperature phase and the 
ground state is striped for arbitrary small dipolar interaction. 
Within the present approach we have gone beyond the usual 
mean field approximation for the stripe phase, showing that 
an isotropic- nematic phase transition takes place at a higher 
temperature than the mean field isotropic-stripe transition. 

Comparing with experimental evidence, as discussed in the 
Introduction, a nematic phase without positional order is prob- 
ably present in ultrathin ferromagnetic films, although it has 
not been characterized already^ 3 -. Our results indicate that the 
nematic phase is more robust for higher values of the ratio 
between the exchange and dipolar interactions, J/g, as is the 
case in experiments where the intensity of the dipolar term is 
two or three orders of magnitude weaker than the exchange 
interaction. 

Comparison with numerical results from computer simula- 
tions is still difficult for several reasons. The first one is that 
our calculation is of mean field character and then our results 
can only be qualitatively correct. A second problem is that up 
to now there have been very few attempts to quantify and char- 
acterize orientational order in systems of the type considered 
in this work, even in computer simulation s 30,32 " ^^ To our 
knowledge, the most detailed simulational study of the phase 
transitions in the DFIF presented so far is the work by Can- 
nas et. ali 2 ^, where evidence was shown of an intermediate 
nematic phase, between a paramagnetic and a stripe-crystal 
phase. It that reference, small values of the ratio J/g were 
considered. In fact, the nematic phase was reported for a ra- 
tio J/g = 2, which was the larger value considered in that 
work. For that case, the nematic phase was observed in a nar- 
row temperature interval, in qualitative agreement with our 
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phase diagram of figure [3] Note also that for J/g = 1 no 
evidence of nematic phase was reported in that work, again 
in agreement with our results, which point to the absence of 
nematic phase for J < 1. It it still a major challenge for 
computer simulations to attain ratios in the experimental range 
J/g ~ 10 2 — 10 3 . To our knowledge, the largest values at- 
tained in simulations have been around J/g ~ 10, which can 
only yield a narrow nematic phas o 34,41 . 

As already said, the present approach is valid for arbitrary 
microscopic interactions, which enter only in the computation 
of the final self-consistent equations. Because of its gener- 
ality, we expect it can be useful for facing some important 
yet unsolved problems regarding the behavior of systems with 
competing interactions. One of these problems is the role 
of the range of interactions in producing pure orientational 
phases. The necessity of long range interactions is frequently 
invoked, but the actual influence of the relative range of the 
competing interactions is still an open problem. 



Some points have to be addressed in order to turn the 
approach quantitative. The computation of the nematic or- 
der parameter is equivalent to the computation of nearest- 
neighbors correlation functions. Then, better approximations 
for the computation of correlation functions will turn the re- 
sults quantitatively reliable. Another point is the formal prob- 
lem with the non-positive character of the quadratic form, 
necessary for applying the Hubbard-Stratonovich transforma- 
tion. Some previous works have already addressed this ques- 
tio n 11 ' 42 , which will be a subject of future work. 

Although we have not tried to get quantitative results, we 
hope the present approach will be useful to describe in full 
the phase transitions in two dimensional systems which show 
microphase separation and nematic-like orientational phases 
originated from competing interactions. 

This work was partially supported by CNPq and FAPERJ 
(Brazil). We wish to thank Sergio Cannas and Alessandro 
Vindigni for a critical reading of the manuscript. 
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